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I. INTRODUCTION 



The relevance of top quark physics at LHC is nowadays firmly established. Exhaustive 
detailed descriptions of the experimental strategies that will be adopted to determine with 
the greatest obtainable precision the observable properties of the measurable processes 
are available in the literature [1]. Within the Standard Model (SM), one sees from a 
careful reading of Rcf. [1] that top production will improve our knowledge of fundamental 
features of the model that are still relatively poorly determined, like the precise value of 
the top mass (aimed LHC accuracy of 1 GeV) or that of the CKM Vtt, couphng (aimed 5 
% precision). For physics beyond the SM, top production offers a variety of possible tests 
that arc also fully listed in Ref. [1]. Generally speaking, these assume the possibility of 
detecting small effects, typically due to virtual exchanges of particles not existing in the 
SM. In order to be able to identify unambiguously such terms, a suitable control of all 
possible (theoretical and experimental) uncertainties is therefore essential. This requires 
dedicated analyses and accurate calculations from all the involved research groups. 

On the theoretical side, a widespread hope exists that LHC will produce a certain 
amount of supersymmetric particles. Assuming that this will be the case, top quark 
physics might then be able to provide consistency tests of the candidate supersymmetric 
model via identification of small supersymmetric virtual effects, provided that the latter 
were sufficiently larger than the intrinsic overall uncertainty that would affect the mea- 
sured process. In the specific framework of the Minimal Supersymmetric Standard Model, 
the calculation of the virtual electroweak effects has been actually performed at the per- 
turbative one-loop level for arbitrary initial partons c.m.energy values. In a first paper 
[2] the production of unpolarized top quarks has been considered. In general, the virtual 
electroweak SUSY effects can become appreciably large, reaching, roughly, the ten per- 
cent size, depending on the values of the several (six, in the approach of Ref. [2]) relevant 
theoretical parameters of the Model. To appreciate and identify these supersymmetric 
contributions, the theoretical uncertainty of the production process, essentially of QCD 
origin, should be therefore controlled and maintained below the 10 % level. This drastic 
operation might be completely avoided if the measured quantity were the QCD-free final 
top longitudinal polarization asymmetry, already considered for LHC measurements in 
a previous paper [3]. In particular, from Rcf. [3] one realizes the existence of a strong 
dependence of the differential asymmetry at LHC on tan /?, for given values of other su- 
persymmetric parameters. Typically, the magnitude of the asymmetry lies in the few 
percent range. To appreciate it, experimental measurements at this precision level would 
be therefore requested. 

The theoretical analyses of Refs. [2], [3] were performed for arbitrary initial parton 
cm. energy values. Under these conditions, all the parameters of the model enter, in 
principle, in the theoretical expressions. To reduce the number of parameters that effec- 
tively determine a prediction would be a welcome feature for any accurate test. In the 
specific case of the MSSM, it has been recently shown [4], that a great simplification can 
be achieved when the final top-antitop pair is produced at a cm. energy ^/s sufficiently 
larger than the heaviest mass M of the SUSY particles that contribute via virtual ex- 
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changes. Assuming a moderately light SUSY scenario, e.g. M ~ 350 GeV, imposes a 
qualitative range ^/s ~ 1 TeV (i.e. s/M"^ ^ 10). Under these conditions, an asymptotic 
logarithmic expansion at one loop of so called Sudakov type can be used to provide a 
satisfactory description of the differential cross sections [5] . To the next-to leading hnear 
logarithmic order accuracy, the only supersymmetric parameter that enters the expansion 
is tan/5 = ^2/^1, via supersymmetric effects of Yukawa kind that only affect the linear 
logarithmic term [5]. All the remaining SUSY parameters are shifted in the next terms 
of the expansion, and in this preliminary analysis we shall neglect them on the basis of 
previous studies performed for LC physics [6], devoting a more rigorous discussion of 
their possible effects at LHC in the final conclusions. 

The aim of this paper is that of investigating which useful information on the MSSM 
would be derivable, in the moderately light SUSY scenario that we described (or in an 
essentially similar one), by accurate measurements of the large (about 1 TeV) invariant 
mass distributions of the unpolarized top-antitop production cross section and of the 
aforementioned final top longitudinal polarization asymmetry. With this aim, we shall 
organize the paper in the following way. In Section 2, we shall quickly exhibit the relevant 
asymptotic expansions for polarized and unpolarized top production, working in the cm. 
frame of the elementary initial partons (essentially, gluon-gluon frame), and show the 
special and simple dependence on tan [3 that appears, discussing qualitatively its possible 
virtues, particularly in the large (tan (3 > 10) values region. In the following Section 3 
the expected theoretical and experimental uncertainties will be discussed, starting from 
the very initial request of transforming the initial partons cm. energy y^-dependent 
predictions into measurable {i.e. top-antitop invariant mass-dependent) expressions and 
making a conservative estimate of the reasonable experimental accuracies obtainable in 
the measurements. Section 4 will show which possible information on the MSSM would 
be derivable under the previous assumptions, providing a final discussion with a number 
of possibly general conclusions. 



II. DIFFERENTIAL CROSS SECTIONS AT PARTONIC LEVEL 

In this Section, we shall perform the preliminary logarithmic expansion of the relevant 
differential cross sections at the partonic level, for production of a polarized top-antitop 
pair, in the cm. system of the initial parton pair. As we anticipated, we shall assume 
a preliminary production of supersymmetric particles, and a moderately light SUSY sce- 
nario where all supersymmetric masses lie below, say, 350 GcV. Under this assumption, 
a reasonable choice for the validity of an asymptotic expansion in the cm. energy ^/s 
seems to us to be ~ 1 TeV, so that, denoting with M the mass of the heaviest su- 
persymmetric particle that provides virtual contributions to the process, s/M'^ ~ 10 (for 
reasonably different values of M, this qualitative request fixes the corresponding value of 

The choice of \fs generates useful simplifications in the theoretical description of the 
process, that have already been stressed in a previous paper [4]. In fact, for production 
at LHC, one already knows that the largely dominant contributions come from the initial 
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gluon-gluon state. At the Born level, the corresponding scattering amplitude is given 
by the sum of three s, t, u channel diagrams. When ^/s is much larger than the top 
mass, only the two t, u channels contribution for opposite helicity gluons survives, while 
the sum of the s,t,u channel diagrams with equal helicity gluons vanishes as ml/s . An 
analogous simplification appears for the final top-antitop pair since the final top helicity 
A = At is opposite to the antitop one Af, neglecting m^/s terms; note that this corre- 
sponds to chirality conservation along the fermionic line since, for the antitop, chirality is 
opposite to helicity. As a consequence of these facts, we shall concentrate our attention, 
at the electroweak one-loop order that we shall consider, on the MSSM corrections to 
the gluon-gluon t,u channel Born diagrams for opposite top-antitop helicities (actually, 
in Ref. [4] a detailed calculation of the one-loop corrections to the initial qq parton state 
was also performed and shown, and the numerical overall calculation confirmed the large 
dominance of the two gluon channel) . 

In the chosen picture, the scattering amplitude for the process is, as we said, only 
given by the sum of the t,u channel terms 

A^-"'* = - ^^m^{Yel){Y{k' -P'^'),){Y^M^ (2.1) 



„2 A-^'A* 

A^-"''' = - ^-^m^{l'^i){Y{k^ - f')pW^>m (2.2) 

whose contribution to the helicity amplitude F^'"'"(/x, — /x. A, —A), where fi and A denote 
the gluon and top helicities, is : 

^Ftr^r., . >x 9,A*A'\ 2Acos^^ + u ^ „ , A-^A\ 2A cos^ + u „ , „, 

F-»fc-^,A,-A) )^-^sm9 + ,J(— )^^-/.nO (2.3) 

where Qs is the QCD coupling constant, i and j refer to the two gluon color states. A* and 
A-' are color matrices, and e{, the polarization vectors, is the angle between the top 
quark and the gluon. 

The corresponding expression of the differential cross sections for fixed top hehcity (L, 
R denote the chirality of both top and antitop) are then : 

dcosO dcosO 8s ^ 3ut As^ ^ ^ ' ' 

The electroweak corrections turn out to be extremely simple at the one-loop level. In 
the adopted logarithmic expansion of Sudakov type, and using the conventional definitions 
of the various terms that one can find e.g. in Ref. [4], only universal gauge and Yukawa 
terms appear. Graphically, they are related to the diagrams shown in Fig.l. Briefly, their 
contribution to the differential cross sections can be expressed by the following simple 
equations : 
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d cos 6 d cos 9 ^ 

_ Tial + _ 3(^^ + t^) 



'm\{\ + cot^ /3) + m^(l + tan^ log (2.5) 



'-f (^^ ^ i^i^) _ rfa^-"(^^ ^ td^^^ ^ 24) 



(i cos Q d cos ^ 

[1 + 



8s ' 3Ki 4s2 

using the universal coefficients c^^ defined in [4,5]. 

Eqs.2.5,2.6 are only intermediate expressions that allow to derive the two quantities 
that we shall consider in what follows. The first one is the (usual) unpolarized top-antitop 
cross section, that in the chosen scenario reads : 

dal'"'^ def da' '°°i\gg ^ + tRh) _ da^'^^^igg ^ + tRtn) ^ 
dcose dcos0 dcosO ^ ^(^l^Cr) 



The second quantity that we shall consider in this paper will be called the final top 
longitudinal polarization asymmetry and denoted by us with the symbol : 

^ da^ '°'^{gg ^ tLh) _ da^^^^^igg-^t^ 

^ ^ dcos^ _ dcos^ _ ^' 

. da' '""Pjgg ^ tLh) da' '""^{gg ^ tnh) . ,^ . 

^ dcos9 dcosd ^ ^ ' ' 

At Born level, it is vanishing. Its theoretical expression at first order in the chosen 
configuration is : 

a, ~ c« - c« - ~ ^^"^^ \2 loe ^ - loe^ ^1 

" STT^M^ ^^^g ]i^]K(l + tan^ /3) - + cot^ Z^)] (2-9) 

It should be mentioned that we have considered a quantity already defined in previous 
papers, in particular for LHC by Kao and Wackeroth [3]. Our definition differs by their 
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definition only by a sign. The reason why wc consider this asymmetry a potentially 
interesting observable has been already indicated in Ref. [3], since at this perturbative 
order at is completely free of QCD effects (and related theoretical uncertainties), and we 
shall come back to this point in the next Section. 

In Eqs. 2.5-2.7,2.9 we have given Sudakov expressions for the one-loop corrections 
in the MSSM case. The corresponding expressions in the SM case can be obtained in a 
straightforward way by the replacement [5]: 2 log — log^ —>■ 3 log — log^, 2m|(l-|-cot^ /?) 
m| and 2m^(l + tan^ (3) ml. 

An important point that we should now clarify is in fact our treatment of the QCD 
corrections to the considered process. Generally speaking they are of two natures, those 
of SM origin and the extra ones of SUSY nature. Since the aim of this paper is the 
identification of genuine supersymmetric effects, the SM QCD effects will be considered 
by us as "partially" known terms, whose theoretical uncertainty will be discussed in 
Section 3. The treatment of the SUSY QCD corrections will be different. For these 
terms, another benefit of the large choice is the disappearance (at the one loop order) 
of the correction to the s-channel diagram, containing the SUSY contributions to the 
gluon self-energy. The remaining term arises from the t, u channel vertex effects with 
gluino exchanges. In our scenario which assumes a gluino mass smaller than about 350 
GeV, its Sudakov logarithmic (linear) contribution to the scattering amplitude is given 
by the following expression: 



-F^-"(/.,-A.,A,-A)^log— (2.10) 

Its effect in the unpolarized differential cross section can be easily computed, and adds 
a contribution that turns out to be, numerically, smaller than that of electroweak origin 
in our scenario. As stressed in Ref. [4] this surviving contribution has the same negative 
sign as that of the electroweak one, which enhances the SUSY effect, and a numerical 
size of a few percent in the considered ^/s region, which makes its one loop expression 
reasonably safe. On the contrary, the suppressed s-channel contribution would have been 
positive, thus decreasing the overall SUSY effect. 

Until now, our analysis has been limited to the description of the elementary partonic 
process, in particular of its electroweak properties. Clearly, this is only a first step towards 
a reahstic prediction for observable quantities, that fully takes into account the effective 
nature of the initial and of the final states. This will be tentatively done in the next 
two Sections. But before entering this concrete stage, we will illustrate with a few extra 
analyses at the partonic level the reasons why we believe that suitable measurements of 
top-antitop production in our proposed scenario might be interesting. With this aim, we 
shall retain for the moment as variable quantity of the process the initial gluon pair cm. 
energy ^/s (whose relationship with the observable invariant mass Mu will be examined 
in Section 3) and assume that an experimental measurement can be performed for two 
proposed distributions. The first one is : 

da{PP ^tt + ...) 1 fcosema. da.j^tt, 

ds = ^ icos.^ f § ^--^^'^^^^^ dW^^^) ] ^2.11) 
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where t = ^, and (ij) represent all initial qq pairs with q — u,d, s, c, b and the initial gg 
pairs, with the corresponding luminosities 

L,,(r, cos^) = — /_ dy [ i{x)3{-) + 3{x)i{-) ] (2.12) 

where S is the total pp cm. energy, and i{x) the distributions of the parton i inside 
the proton with a momentum fraction, x = e^, related to the rapidity y of the tt 
system [7]. The parton distribution functions are the 2003 NLO MRST set available 
on [8] . The limits of integrations for y can be written 

Vmax = max{0, min{F - ^ log x, ^ + ^ ^og x, - log( V^)}} 

ymin — Umax 

(2.13) 

where the maximal rapidity is y = 2, the quantity x is related to the scattering angle in 
the tt cm. 

1 + cos^ . ,N 

1 — cos 9 

and 



cos (7mm, max — Tu ~ {■^■^'^) 



expressed in terms of the chosen value for pT,min — 10 GeV. 
In an analogous way, we define : 



A,is) = - ^]/[^ + ^] (2.16) 
ds ds ds ds 



with 

— = - / d cos e [ Lij{T, cos 9) — „ s 2.17 

ds b J cos dmiu T; dcosU 

Wc have drawn in Figs. 2 and 3 the values of the two distributions versus the cm 
energy ^/s in the energy range that would be relevant in our scenario, from roughly 700 
GeV to roughly 1 TeV, for three significant values of tan/3 i.e. tan/3 = 1, 10,50. These 
have been obtained integrating our expressions of the two gluon process written in this 
Section and adding to them the extra depressed initial quark-squark terms, given in Ref. 
[4] and not rewritten here for simplicity purposes. As one sees, At{s) changes sign in the 
range. It should be stressed that our values of this asymmetry are in excellent agreement 
with the previous results of Ref. [3] (Fig.3 of that paper) if the comparison is limited to 
the set of parameters of that Reference that would correspond to our scenario (i.e. those 
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corresponding to the smallest SUSY masses). One sees also that the electroweak MSSM 
effect could he relatively large. In da/ds it would reach a relative 20-22 % for tan/? ~ 50 
and ^/s = 1 TeV ; in At{s), the effect is smaller, of the six percent size for the same values 
of the parameters. One should keep in mind, though, that it represents in this case, not 
a small correction to a main SM value, but essentially the overall value of the asymmetry 
(which would be at the one percent level in the SM). 

An additional feature of the two considered distributions is that the dependence on 
tan/3 is rather strong. This can be understood by looking at Eqs. (2.7-2.9) and is graphi- 
cally shown in Figs.4-5. One notices in particular an enhancement of the effect for large 
tan/3 values in both observables. In the asymmetry the effect changes sign when tan/3 
moves from the low (~ 1) values to the large (~ 50) ones, a fact that can also be easily 
understood looking at Eq.2.9. 

Encouraged by these promising features, we considered at this point at a qualitative 
level the possibility of deriving confidence limits on tan/3 from a analysis of the two 
observables. With this purpose, we have used the expected event distribution shown in 
Fig.6 and taken from [9]. It has been obtained from the semileptonic decay channel into 
muons only. The integrated luminosity, taking into account the branching ratio, is of 
about 5 fb~^, corresponding to half a year of data taking at low luminosity (lO^^cm^s"^). 
This has provided us with a realistic statistical experimental error and for this analysis 
we have not considered the other sources of uncertainties, that will be discussed in the 
next Section 3. Note that the event distribution that we have used is given as a function 
of the top-antitop invariant mass Mu. In this qualitative analysis, we have assumed the 
equality Mu = ^/s, whose validity will also be discussed in the next Section 3. 

In agreement with the spirit of this preliminary paper, we shall retain, we repeat, the 
two (leading and next to leading) terms of the logarithmic expansion in the analysis. A 
discussion of the possible effects of next to next to leading terms will be given in the final 
conclusions. 

The results of our fits to tan/3 are shown in Fig. 7, in the two cases of (a) measure- 
ment of only dau/ds and (b) measurement of both dau/ds and At. As one sees from these 
figures, it would be possible to obtain from our analysis the following main information: 

1. A measurement of only dau/ds would provide information in two separate ranges, 
1 < tan /3 < 3 and 20 < tan /3. This can be understood as follows. In our approxi- 
mation, the unpolarized cross section depends on tan /3 through the precise Yukawa 
correction which is a definite combination of tan^ /3 and cot^ (3. Each value of the 
correction is obtained from a couple of mirror values of tan /3. A analysis is thus 
plagued by this irreducible uncertainty and will always show a false minimum at 
the mirror value of the true tan/3. However, if 1 ^ tan/3 ^ 3 or 20 ^ tan/3, there is 
a Ax^ = 1 confidence region centered around the exact value of tan /3 and we have 
drawn in the Fig. 7 the corresponding boundaries. If, on the other hand, tan/3 falls 
in the dead region 3 < tan /3 ^ 20, it is practically impossible to tell between the 
true and false (mirror) minima of x^ and we have not drawn any line. 

2. The additional information obtained by the measurement of At is crucial to eliminate 
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the above mentioned dead region 3 < tan /3 < 20. 

3. The more interesting information is, in our opinion, that relative to the large tan /S 
region, roughly tan/3 ^ 20. In fact, in the low tan/3 ~ 1 region there should be 
alternative accurate determinations provided by other measurements [10]. On the 
contrary, for the large tan/3 values, no accurate alternative measurements at the 
LHC time seem to be, to our knowledge, available. 

4. In the large tan/3 region, while the measurements of dajj/ds alone seems to provide 
already a satisfactory result (with a relative error below 20 % for tan /3 larger than 
about 30), the addition of At would lead to an extremely accurate determination 
reaching a few percent value for tan /3 ~ 50. 

This quahtative analysis is now completed. We feel that the features that have emerged 
justify the effort of moving to the detailed consideration of the related physical process, 
and of examining whether the nice properties that we underlined at qualitative level 
survive in the realistic experimental measurements. This will be done starting in the 
forthcoming Section 3. 

III. EXPERIMENTAL AND THEORETICAL UNCERTAINTIES 

The discussion of virtual supersymmetric effects has led to corrections to the lowest 
level of the partonic cross section for top-anti-top production as function of the partonic 
centre-of-mass y^. Notably in the high region of ^/s, the corrections can be as large as 20 
to 30% and depend on the exact value of the supersymmetric parameters. In this section 
we concentrate on the experimental issues involved in order to measure these corrections. 

Although it can be inferred from Fig. 2 that the shape of the distribution of M^t as 
function of is distorted by the supersymmetric effects, most of the sensitivity to these 
new effects originate in the difference of the ti cross section above y/s of approximately 600 
GeV. The determination of absolute cross sections at hadron machines is however highly 
non-trivial. The estimated overall predicted theoretical uncertainty of tt production is of 
the order of 12% [1] , the error coming from remaining scale dependence of the NLO QCD 
calculations. 

Experimentally there are two main concerns to determine the ti cross section as func- 
tion of \fs. First, one can only measure the final state top pairs, and has therefore 
experimental access to the invariant mass Mu only. NLO QCD effects (final state gluon 
radiation, virtual effects) spoil the equivalence of Mu with ^/s. Second, detector effects 
like efficiency uncertainties, jet resolutions and mis-calibrations migrate the true Mtt dis- 
tribution. We discuss now these two effects. 

A. Higher order QCD effects 

Higher order QCD effects perturb the determination of the cross section as function 
of in two ways. The tt cross section increases from 590 pb to 830 ± 100 pb, calculated 
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in LO and NLO respectively, at the LHC. Besides the normahsation, also the shape of 
the Mft distribution gets distorted by NLO effects, due to real and virtual effects. 

The effects of the NLO QCD calculations have been investigated using the 
MC@NLO [12] program, that incorporates a full NLO treatment in the Herwig MC gen- 
erator. We have used this program to generate the Mu distribution both in LO and in 
NLO, using identical parton shower, parameters settings etc. The effect of the NLO cal- 
culations on the shape of the Mu distribution is shown in Fig. 8, where the ratio NLO/LO 
is given as a function of M^f at NLO. The value of M^t is obtained at the parton level, 
as the invariant mass of the top and anti-top quark, after both ISR and FSR. The LO 
and NLO total cross sections were normalised to each other. Deviations from unity in 
this figure are entirely due to differences in shape of Mu- As one sees from the figure, the 
relative difference between and Mu remains bounded (below roughly 5%) when ^/s 
varies between ~ 700 GeV and 1 TeV, which is the chosen energy range of this paper. 
For larger ^/s values the difference raises and can reach a 10 % limit when approaches 
what we consider a reahstic experimental limit (see Fig. 6) of the search, i.e. ^/s ~ 1.3 
TeV. 

B. Experimental systematic uncertainties 

Top events are triggered and selected best when they decay semileptonically, i.e. where 
one top decays to a b-jet plus a lepton and a neutrino, the other top decays to a b-jet and 
two light quark initiated jets. The lepton is ideal to 'trigger' the event, and the two b-jets 
are clear signatures for top production. The kinematics of the neutrino can be recovered as 
the missing transverse energy of the event. Ambiguities in the longitudinal direction can 
be dealt with by consistency requirements on the event topology (for example requiring 
two tops to have the same mass). We generated 1.1x10^ ti events, using the Pythia 
Monte Carlo [13] and processing them through the ATLAS detector fast simulation [14] . 
The number of events correspond to about 5 fb~^ of collected data. The selection and 
reconstruction criteria applied are described in Ref. [9], resulting in an efficiency of about 
1.5%. This overall efficiency results from the stringent cuts applied to select the events, 
including the b-tagging. These requirements were applied to ensure a high purity of our ti 
sample, and the surviving events contain a negligible amount of background not originated 
from ttbar production. A first estimate of the experimental errors involved is obtained 
by considering two of the main sources of systematic uncertainties in the determination 
of the invariant mass Mtf, the jet energy scale uncertainty and the uncertainties of jet 
development due to initial and final state showering. To evaluate the effect of an absolute 
jet energy scale uncertainty, a 5% miscalibration coefficient was applied to the jet energies. 
This produces a bin-by-bin distortion of the distribution smaller than 20%. This is 
certainly an overestimate of the possible error, since it has been shown that in ATLAS an 
accurate absolute energy calibration of light quarks and b-jets can be extracted from Z+jet 
events, with an expected precision of about 1% [1]. An in-situ calibration of light jets in 
which both the absolute energy and direction cahbration are extracted from the W — > jj 
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channel itself is possible as well. For this purpose a cleaner sample of W candidates can 
be selected from the ti events. 

Uncertainties in the initial state radiation from the incoming partons (ISR) and final 
state radiation from the top decay products (FSR) affect the precision of the M^f measure- 
ment. To estimate this effect, as suggested in [1], the Mu distribution obtained with the 
standard Pythia Monte Carlo has been compared with the same distribution determined 
with ISR switched off. The same approach was employed for FSR. The level of knowledge 
of ISR and FSR is of the order of 10% so the systematic uncertainty on each bin of the 
tt mass was taken to be of 20% of the corresponding difference in the number of events 
obtained comparing the standard mass distribution with the distribution having switched 
off ISR (or FSR). This results as well in an error that is smaller than 20 %. 

Experimentally, the luminosity introduces an ultimate uncertainty of the order of 
5%; at the startup period of LHC this uncertainty will be much larger. In conclusion, 
from a rather conservative evaluation of the most relevant, theoretical and experimental, 
uncertainties, we are led to the conclusion that an overall error of approximately 20 %- 
25 % in our energy range appears realistically achieavable from a preliminary estimate. 
This does not exclude the possibility, that appears to us to be strongly motivated, that 
further theoretical and experimental efforts might reduce this value to a final limit of, say, 
15 or 10 % size. In the final Section of this paper, we shall show what might be the level 
of theoretical information obtainable under these assumptions. 

To conclude this Section, we want to discuss briefly the possibility that we mentioned 
of using as realistic experimental input the quantity At, the final top longitudinal polar- 
ization asymmetry, defined in our Eq. (2.16), whose potential virtues were qualitatively 
investigated in Section 2. As we said in that Section, this quantity would exhibit the 
remarkable property of being essentially free of QCD effects (and, more important, of the 
related theoretical uncertainties). In our opinion, this would motivate a rigorous experi- 
mental analysis, like the one that we have shown for the unpolarized cross section. This 
analysis is not, at the moment, available. To provide a motivation for it, we shall show 
in the next Section the additional information that would be achievable if both measure- 
ments of au and Af were performed. With this purpose, we shall simulate for A^ three 
different overall errors of the same relative size (i.e. 20, 15, 10 %) that we assumed for 
au- Although we realize that these values might be unduly optimistic (or, possibly and 
hopefully, pessimistic), we believe that the improvements of information that would be 
obtained should justify a dedicated experimental effort in that direction. 

Section 3 is thus concluded. The following short Section 4 will be devoted to a review 
of the potential information that would be derivable from the results that we have listed, 
added to the formulae that we have derived in Section 2. 

IV. INFORMATION ON tan (3 FROM A REALISTIC ANALYSIS OF THE 
LARGE INVARIANT MASS DEPENDENCE OF TOP OBSERVABLES 

In this final Section we shall show the results of the generalization of the analysis 
performed in Section 2, under the assumption of a purely statistical error, to the more 
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realistic cases of an additional overall systematic (theoretical + experimental) relative 
error on the unpolarized cross section. The results of the fits are shown in Fig. 9. In the 
spirit of the last remarks of Section 3, we have also computed in Fig. 9 the results of three 

fits where, in addition to the unpolarized cross section, the longitudinal polarization 
asymmetry At has also been included, assuming for it a relative error equal to that of the 
unpolarized cross section. 

In more details, the left side of Fig. (9) shows the relative accuracy on tan /? that can 
be derived by exploiting both au and A-t and adding to the statistical error a systematic 
error constant in energy with three possible values, 10%, 15% and 20%. The right side of 
the Figure shows the same analysis, but without exploiting Af. From a comparison of the 
two plots, one sees that the contribution from At is definitely not negligible especially at 
larger values of the systematic error. For instance, the determination of tan f3 using both 
au and At with a systematic 20 % error gives an accuracy roughly equal to that obtained 
using only au with a systematic 15% error. 

Apart from the specific role of At, we stress that the results shown in Fig. (9) are, in 
our opinion, encouraging. Indeed, even in the worst case, i.e. discarding At and using 
only au at 20%, we find an accuracy on tan/9 that is around 35% at tan/? ~ 50. We 
emphasize again that this is due to the known special role of tan/3 appearing directly 
in the couplings and allowing a wide range of enhancement in the radiative corrections. 
Indeed, even with an error of 20% on the experimental measurements, the radiative effects 
at tan /3 > 50 can be large enough to be visible and able to constrain the possible values 
of tan (3 in the fit. 

We also stress, to conclude this Section, that the hypothesis of a systematic error of 
20 % seems to us indeed pessimistic since, for instance, an overall 5% is the declared goal 
required for a determination of mt with a precision of 1 GeV ( [1], page 428). 

V. CONCLUDING REMARKS 

As a final comment to our work, we would like to make the following statement. The 
main purpose of our paper was to show that, at LHC, precision tests of the Minimal 
Supersymmetric Standard Model would actually be possible, even if the experimental 
conditions will never be at the level of precision of LEPl, where the Standard Model was 
tested at the "less than one percent" level. The qualitative level of "optimal" accuracy 
would typically be here at the "ten percent" level, perhaps a little worse but possibly a 
little better. Under these conditions, we have shown that, even at an accuracy level of 
(only) twenty percent, the production of top-antitop pairs in a region of relatively large 
(but experimentally realistic) invariant mass could provide a non trivial confirmation (or, 
perhaps, a prediction) for the value of one of the most important parameters of the model, 
tan/5. The technical reasons of this possibility are that (a) this process contains contri- 
butions of Yukawa kind from vertices where the heavy top mass enters quadratically, 
enhanced by a large Sudakov logarithm and by quadratic factors of tan P and (b) that 
tan P is actually the only SUSY parameter that enters the logarithmic expansion. The 
combination of these enhancements might in fact lead to an observable effect, in particular 
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from a proper fit to the invariant mass distribution. Thus, from the accurate measure- 
ments of the production of the heaviest quark pair producible in the reaction, important 
tests of the model (or even predictions, in case tan (3 were still poorly known) might be 
derived. We like to remark that an analogy would exist with the analogous situation met 
at LEPl, where, from accurate measurements of the producion of bb pairs, fundamental 
tests and "predictions" for the value of the top mass were derived. 

One should also (honestly) repeat, to conclude, that our (preliminary) analysis was 
performed assuming special conditions both for the model (a "moderately' light SUSY 
scenario) and for the logarithmic Sudakov expansion (next-to -next to leading order). In 
a less favourable case, though, our analysis could be repeated with the proper values of 
all the SUSY parameters and a complete one-loop description, which would not benefit 
from the simple Sudakov expansion. This more complete analysis in in fact already being 
carried on in [15]. 
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FIG. 1. Diagrams for electroweak Sudakov logarithmic corrections to gg — > tt. In (a),(b) the 
virtual gauge boson can be a photon, a Z or a W^. In (c) the dashed line represents a Higgs 
or a Goldstone particle {H^'^, h^, A^, G^'^). In (d) the dashed lines represent a sbottom or a 
stop and the solid line inside the triangle represents a chargino or a neutralino. 
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FIG. 2. Relative one loop corrections to the unpolarized cross section dau/ds. More pre- 
cisely, wc show the difference d(j]j^""^/ds - dcj^°"'/ds divided by da§"''''/ds. The plot shows the 
effect in the Standard Model and in the MSSM with three reference values of tan (3 = 2, 10, 50. 
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FIG. 3. One loop value of the top quaxk polarization asymmetry At. The plot shows the 
effect in the Standard Model and in the MSSM with three reference values of tan (3 = 2, 10, 50. 
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FIG. 4. Relative one loop corrections to the unpolarized cross section dau/ds. As in Fig. 2, 
we show the ratio {da]j^°'^ / ds - du^'^'' / ds) / da^""'' / ds . The plot shows the effect in the MSSM 
as a function of tan (3 at two reference values of ^/s = 700, 1000 GeV. 
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FIG. 5. One loop expression of the top quark polarization asymmetry Af. The plot shows 
the effect in the MSSM as a function of tan /? at two reference values of s/s = 700, 1000 GeV. 
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FIG. 6. Event distribution as function of M^j in the Standard Model. The events are obtained 
from scmilcptonic top decay corresponding to 5 fb~^ of collected data (half a year at luminosity 
lO^^cm-^s-i). 
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FIG. 7. Results from the analysis of the tan/3 dependence of the two considered distri- 
butions. For each hypothetical true value of tan (5, the figures show the Ax^ = 1 boundaries 

tan/3+ A-t and the corresponding relative error defined as ^(A+ — A—)/ tan /3. The dashed and 
solid lines arc obtained by exploiting the unpolarized cross section alone or the combination of 
it with the polarization asymmetry. 
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FIG. 8. Ratio of normalised NLO and LO predictions of the distribution of Mu, i.e. 
\Mtt/J2Mtt]/[y/s/J2 \/^) determined from the MC@NLO and Herwig programs. This ratio 
is only sensitive to the difference in shape of the two predictions. 
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FIG. 9. This figure shows the same analysis described in Fig. (7). However, we now take 
into account an energy independent fixed systematic error of 10, 15 or 20 %. The curves on 
the left are obtained by including in the analysis both the unpolarized cross section and the 
polarization asymmetry At. The curves on the right do not use Af. 
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